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ABSTRACT 

We analysed high-resolution UVES spectra of six stars belonging to the subgiant branch of o) Centauri, and derived abundance ratios 
of 19 chemical elements (namely Al, Ba, C, Ca, Co, Cr, Cu, Fe, La, Mg, Mn, N, Na, Ni, Sc, Si, Sr, Ti, and Y). A comparison with 
previous abundance determinations for red giants provided remarkable agreement and allowed us to identify the sub-populations to 
which our targets belong. We found that three targets belong to a low-metallicity population at [Fe/H]^-2.0 dex, [a/Fe]^+0.4 dex and 
[s/Fe]^0 dex. Stars with similar characteristics were found in small amounts by past surveys of red giants. We discuss the possibility 
that they belong to a separate sub-population that we name VMP (very metal-poor, at most 5% of the total cluster population), which 
- in the self-enrichment hypothesis - is the best-candidate first stellar generation in lj Cen. Two of the remaining targets belong to 
the dominant metal-poor population (MP) at [Fe/H]^~1.7 dex, and the last one to the metal-intermediate (Mint) one at [Fe/H]=:- 
1.2 dex. The existence of the newly defined VMP population could help to understand some puzzling results based on low-resolution 
spectroscopy (Sollima et al., Villanova et al.) in their age differences determinations, because the metallicity resolution of these studies 
was probably not enough to detect the VMP population. The VMP could also correspond to some of the additional substructures of 
the subgiant-branch region found in the latest HST photometry (Bellini et al.). After trying to correlate chemical abundances with 
substructures in the subgiant branch of a) Cen, we found that the age difference between the VMP and MP populations should be 
small (0±2 Gyr), while the difference between the MP and Mint populations could be slightly larger (2+2 Gyr). 

Key words, stars: abundances - stars: main sequence -globular clusters: individual {lu Centauri); NGC 6397 



1. Introduction 

Much has been written on co Cen, which is probably the mos t 
studied cluster in the Milky Way (Ivan Leeuwen et all l2002h . 
From the pioneering studies in the 60s to the latest high-quality 
data and models, more and more details of its multiple stellar 
populations have come to light, but the picture did not become as 
clear as expected. Excellent photometries and astrometric cata- 
log ues have r e cently been produced by both ground-based (such 



l2005bl [Stanford et al.L 120061: IVillanova etaLl l2007l) resolution. 
Nevertheless, several puzzles still await solutions. 

One of the open problems concerns the most sensitive re- 
gion of the colour magnitude diagram (hereafter CMD) to 
age differences: the subgiant branch (SGB). A large num- 
ber of photometric and low-resolution sp ectroscopic studies 
(iHughes & WaUersteinLl2000HHilker & Rich tler. 200Q^anc^ 



tiy 

LeeetalL Il999t iPancino et al.L l2000t Ivan Leeuwen et al 



2000t iHilker & Richtleri l2000t iHughes & WallersteinL I200C : 



Sollima et al.l 
20091) and 



I2005al: Icalamida et all 120091: iBellinietal 



Bedin et all 
by high-quality 
high (e.g.. 



space 
12004 



t elescopes (e .g., Ferraro et all 2004 



,2003; 'Hughes et al.', '2004' 'Hilker et al., 2004t iRey et alll2004 
Ferraro et al.^ i2004i : iStanford et al.. 2006i) found age spreads 



iBellini et al.L 120 101) . complemented 
spectr oscopic abundance studi es both at 
Norris & Da Costa. 1995; Smith et aL[ 
Cunha et al., 2002; Pancino, 2003; Johnson et al, 2008! 
Johnson & Pilachowski .2010i: .Marino et al.i, 201(1) and low 
(iNorris et all Il996t ISuntzeff & KraftL ll996ns'ollima et al.L 



200C; 



2009; 



* Based on data collected at the ESO VLT in Chile, with UVES and 
FLAMES under programs 68.D-0332(A) and 079.D-0021. Also based 
on literature data from the ESO WFI, under programs 62.L-0354 and 
63.L-0439, and on data collected from the following online databases: 
NIST, VALD, Kurucz, GEISA. 



ranging from 2 to 6 Gyr, with a few exceptions and puzzles (see 
Section 15.3! for more details). One example of the difficulties 
encountered i n the study of t his com ple x region is posed by th e 
two studies bv lSollima et alJ(l2005bl) and l Villanova et al.l(l2007l) . 
who used the same ACS dataset and low-resolution spectra of 
similar quality, but reached opposite conclusions on the total age 
spread - and age distributio n - of to Cen. Still, even with the best 
ACS photometries (see, e.g. lBellini et al.L 1201 01 it is not easy to 
understand which features of the SGB region correspond to each 
of the populations that are spectroscopically identified on the red 
giant branch (RGB), which are known in great detail thanks to 
the above cited works. This understanding is crucial to solve the 
relative ages problem in ai Cen, and to derive the age-metallicity 
relation, a fundamental ingredient of any model for the forma- 
tion and evolution of this unique stellar system. 
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Table 1. Observing Logs. 
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R.A. (J2000) 
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(B-V)o 


(V-I)o 


R 


(for) 
'-exp 


S/N 


V,- 




Notes 




(deg) 


(deg) 


(mag) 


(mag) 


(mag) 




(sec) 


(@550nm) 


(km s" 


') 




503358 


201.960086 


-47.759426 


17.53 


0.65 


0.65 


45 000 


10800 


50 


231.0 ± 


0.2 


Lower SGB 


503951 


201.870201 


-47.746107 


17.58 


0.65 


0.65 


45 000 


11400 


50 


235.8 ± 


0.1 


Lower SGB 


507109 


201.848173 


-47.692322 


17.31 


0.60 


0.66 


45 000 


11288 


50 


229.7 ± 


0.5 


Upper SGB 


507633 


201.813303 


-47.685166 


17.39 


0.56 


0.66 


45 000 


12000 


45 


226.7 ± 


0.3 


Upper SGB 


512115 


201.846338 


-47.637524 


17.64 


0.65 


0.66 


45 000 


11400 


45 


241.4 ± 


0.1 


Lower SGB 


512938 


201.804070 


-47.630694 


17.32 


0.60 


0.65 


45 000 


10000 


50 


234.6 ± 


0.2 


Upper SGB 




Fig. 1. Location of the programme stars on the CMD of 
u) Cen. The WFI B B-I photometry of stars in CCD#5 (from 
iPancino et al.L l2000h is shown as grey dots. The six UVES tar- 
gets are marked (black filled circles) with their WFI catalogue 
numbers. 

In this paper we present the analysis of a set of UVES 
high-resolution spectra of six SGB stars, selected f rom the wide 
field p hotometric catalogue bv lPancino et al ] (120001) andIP ancind 
(l2003h . A preliminary an alysis of the same dataset was pre- 
sented by iPancinol (l2003h . We describe the spectra reductions 
in Section|2j the abundance analysis in Section[3j and the abun- 
dance ratio results in Section |4] Our main results are discussed 
in detail in Section|5]and are summarized in Section|6] 

2. Observations and data reduction 

We selected our six target s from the WFI B and I photometry 
presented in Pancino et al.l (l2000l) . complemented with V mag- 
nitudes from Pancino! (l2003h . Th e coordinates were obtaine d us- 
ing the astrometric catalogue by Ivan Leeuwen et al.l (l2000h . As 
shown in Figures [1] and |2] the SGB region of u Cen in the ex- 
ternal parts of the cluster shows clear substructures (all our tar- 
gets lie on the WFI CC D#5). This is more clearly seen in other 
literature photometries (iFerraro et al 1 l2004t iBedin et al.L l2004t 



Fig. 2. Location of the programme stars on the area of d Cen. 
Grey dots m a rk sta rs belonging to the WFI photometry by 
iPancino et al.l ( l2000l) . Filled circles mark the postion of the six 
UVES targets. 

Be llini et al l l2010l) .butthev are obtained from space, with the 
HST, in the very centre of o) Cen, where UVES follow-up from 
the ground would be difficult. Three of the programme stars were 
selected towards the upper envelope of the SGB and another 
three towards the lower envelope. 

Echelle spectroscopy was obtained on 18-20 March 2002, 
with UVES at the ESO Kueyen (VLT UT2) telescope, on Cerro 
Paranal, Chile. Table [T]reports the log of the observations, along 
with some basic target information. Each star was observed 
twice and on different nights, both to minimize the cosmic rays 
impact and to identify possible radial velocity shifts. Given the 
faint magnitude of these stars (V^17.5), the spectra were binned 
on chip (2x1), so that the resolution element is covered by ^3 
pixels. A final S/N^50 per pixel was achieved around 550 nm. 

The red spectra (upper and l ower red CCDs ) were reduced 
with the ESO-UVES pipeline dBallester et al.L l2000h . which 
semi-automatically performs bias correction, flat-field correc- 
tion, inter-order background subtraction, optimal extraction with 
cosmic-ray rejection, wavelength calibration (with rebinning), 
and final merging of all overlapping orders. However, since the 
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Table 2. Equivalent widths and atomic data for the lines used in the classical abundance analysis of the program stars. The complete 
version of the table is available at CDS. Here we show a few lines to illustrate its contents. 
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4283.01 


Cal 


1.89 


-0.29 


77.0 


8.9 


1.463 


95.5 


12.2 


1.832 . 


0.0 


0.0 


0.000 


67.0 


5.6 


1.248 


4289.37 


Cal 


1.88 


-0.39 


75.6 
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1.566 


80.2 
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1.378 . 


79.7 


12.7 
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64.5 
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4298.99 
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-0.51 
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0.0 


0.000 
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0.0 


0.000 . 
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0.0 


0.000 


0.0 
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0.000 


4302.53 


Cal 
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0.18 
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0.000 
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0.000 . 
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0.000 
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0.000 


4318.65 


Cal 


1.90 


-0.29 
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8.1 


1.876 
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7.5 


1.747 . 


0.0 
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0.000 


61.1 


3.7 


0.851 


4434.96 


Cal 


1.89 


0.07 


0.0 


0.0 


0.000 


0.0 


0.0 


0.000 . 


0.0 


0.0 


0.000 


81.0 


5.9 


1.227 



S/N ratio is significantly lower on the blue part of the spectra 
(S/N^25 around 450 nm), we decided to manually perform the 
echelle reduction for the blue CCD with the noao.imred.ccdred 
and noao.imred. echelle packages within IRAlfl 

The two one-dimensional and wavelength-calibrated spec- 
tra obtained for each star were normalized by fitting their con- 
tinua with a cubic spline, then corrected for the main telluric 
absorption features (noao.onedspec. telluric in IRAF), using as 
a reference a hot, fast rotating star ( HR 5206/HD 120640) se - 
lected from the Bright Star Catalog dHoffleit & JaschekL Il99lh . 
and observed each night at an airmass not too diff'erent from the 
targets. Finally the spectra were summed - after correcting for 
radial velocity shifts - to produce one single spectrum for each 
of the six observed stars. 

Radial velocities were obtained with DAOSPE(j3 
dStetson & Pancinol l2008h and the procedure described by 
iPancino et al.l ( 2010l) . In short, the laboratory wavelength of 
selected absorption lines (Section 13.11 ) was used to measure 
the observed radial velocity. The heliocentric coiTection was 
computed with IRAF and the telluric H2O and O2 absorption 
bands redward of 580 nm were used to correct for zeropoint 
shifts. All six stars were radial velocity members of a> Cen, con- 
sid ering an av erage V ,.=232 .8 or 233.4 km s~', as determined 
by iMevlan et al. (1995) and lPancino et al.l (l2007h . respectively, 
with a central velocity dispersion of the order of 20 km s ' . The 
resulting velocities and errors are listed in Table [1] 



3. Abundance analysis 

3.1. Equivalent widths and atomic data 

We selected the majority o f our Unes and of t heir atomic data 
from the VALE0 database ( iKupka et all Il999l) . To identify re- 
liable lines, high S/N median spectra of the six UVES targets 
were created and the cleanest, unblended lines of the available 
elements were identified. Only lines that appeared in at least 
three of the six stars were retained in our preliminary selection. 
DAOSPEC was used to measure the equivalenth widths (EW) 
of all the chosen lines. A first-pass abundance analysis was per- 
formed (Section 13.31 ): lines that sho wed systematically higher 
errors and bad Q p arameters (see IStetson & Pancinol 120081: 
iPancino et al.L l2010l for details) and that simultaneously gave 
systematically discrepant abundances were rejected. Finally, all 



' http://iraf.noao.edu/. IRAF is distributed by the National Optical 
Astonomy Observatories, which is operated by the association of 
Universities for Research in Astronomy, Inc., under contract with the 
National Science Foundation. 

^ http://www3.cadc-ccda.hia-iha.nrc-cnrc.gc.ca/community/STETS 
ON/daospec/ ; |http://www.bo.astro.it/~pancino/projects/daospec.html| 

' http://www.astro.uu.se/~vald/ 



WFI 512115 



AEW - 0.2 = 4.6 ml 

^ I I I ^ 

40 



Fig. 3. Comparison between DAOSPEC and IRAF splot mea- 
surements for star WFI 512115. 



lines that gave EW<15 mA, or EW>100 mA (with a few excep- 
tions, see Section [33b were not used to determine abundances. 
The DAOSPEC EW measurements used for the abundance anal- 
ysis are shown in Table |2l along with the formal error ^EW an d 
the quality parameter Q for each line dStetson & Pancinol l2008h . 

A few lines were measured with the help of spectral syn- 
thesis because they were reported to have significant hyper- 
fine structure (HFS ). In particular, we used the atomic data by 
iMartin et al.l (fl988l) for the Mn I lines at 4030, 4033, 4034, 4041, 
and 4055 A; the NIST|^ atomic data for the Ball Unes a t 4934 
5853, 6141, and 6496 A; the atomic data bv lLawler et~ai](l2001l) 
for the Lall U nes at 3988, 40 86, 4123, and 4238 A; and the 
atomic data bv iBielskil (Il975h for the 5105 A Cu I lin^ For 
the CH and CN molecular bands, we used the Kurucz molecular 
lineUstfl but we had to revise the logg f values of C downwards 
by 0.3 dex, similarly to what was reported by Bonifacio et al] 
(ll998h . lLucatello et al.l d2003l) . and lSpite etal] d2005h (see also 
Sectionl33Tl. 

3.2. Atmosplieric parameters and best-model search 

A first guess of the atmospheric parameters was derived from the 
WFI photometry. Dereddened (B-V)o and (V-I)o colours were 
obtained fro m B, V, and I magnitudes adopting E ( B-V)= 0. 1 1 
JLuR I2OO2I) and E(V-I)/E(B-V)=1.30 dDean et all [19781) . and 
are listed in Table[T]along with the V magnitudes. The V-I colour 
was converted from the original V-Ic, based on the Cousins I 
magnitude, to the V-I j, based on the Johnsons I magnitude, with 



http://physics.nist.gov/PhysRefData/ASD/index.html 
We could not use the 5872 A Cul line since it falls into the gap 
between the two UVES red CCDs. 
* http://kurucz.harvard.edu/LINELISTS/LINESMOL/ 
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Table 3. Atmospheric model parameters. 




6550 



Wavelength (A) 



Fig. 4. Example of the H„ profile fit for star WFI 512938. The 
solid line represents the best-fitting synthetic spectrum where all 
parameters are kept fixed as in Table [3] except for the effective 
temperature, T^eff- The inset panel shows the run of the;^^^ of fits 
with different values of Tf//. 



the relations by 'BesselJl ( Il979l) . Effective temperatures (here- 
aft er Teff) and bo lometric corrections (BCv/) were obtained with 
the lAlonso et alj (Il999b calibration. Surface gravities (hereafter 
\ogg) were then derived by means of fundamental relations: 



\ogg, = 0.4(Mv + BCv) + 41og T,//., - 12.61, 



(1) 



where the solar values where assumed i n conformity with the 
lAU recommendations dAndersenl [T999h . i.e., log go = 4.437, 
Te//,o=5770 K and Mio,,Q=4.75. A t ypical mass of 0.8 Mp was 
assumed for the programme stars dBergbusch & VandenBerd, 
20011) . and we used (m-M)i/=14.04±0.11 mag (iBellazzini et al.i. 



20041) . An independent estimate of T^//, which is the most in- 
fluential parameter when determining abundances, was also de- 
rived from the profile fitting of the Hq, line wings. We compute d 
Kurucz atmosphere models with ATLAS9 (lKuruczL[T993ll2005h . 
using gravities and global metallicities close to the photomet- 
ric parameters of each target. We built synthetic spectra with 
a modified version of SYNTH^i exported for the Linux OS 
dSbordone et all l2004t) . An example of a typical H„ profile fit 
is shown in Figure|4] 

We found general agreement between the photometric and 
the H„ Teff estimates, but the scatter was large, with differences 
of up to 300 K in some cases. When averaging estimates from 
B-V and V-I together, the scatter went down (see Table O. We 
decided to rely on the H,, temperatures as a first guess for the 
spectroscopic analysis. First guess log g values were derived 
from H p, temperatures using V magnitudes and the lAlonso et alJ 
(119991) calibration. 

The final atmospheric parameters were then derived through 
the usual "spectroscopic method" (see Section [373l for abundance 
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calculation details): after we had a good number of Fe I and Fe II 
lines (approximately 200 and 10, respectively, depending on the 
star), we re-adjusted the parameters by choosing those that min- 
imized ( i) the slope of [Fe I/H] versus the excitation potential, 
Xex\ (ii) the slope of [Fe I/H] versus the EW of each line|3; (Hi) 
the difference between the average [Fe I/H] and [Fe II/H]; and 
finally by checking that [Fe I/H] did not change significantly 
with wavelength. The adopted atmospheric parameters are listed 
in Table |3] along with the photometric and H„ temperatures, for 
comparison. We found that, on average, spectroscopic Teff es- 
timates based on Fe were lower than H„ ones by 17+52 K, and 
also lower than photometric ones by 83+26 K. 



3.3. Abundance calculations 

For most chemical species we computed abunda nces with th e 
help of the updated version of the original code by'SDite'('l96 



Our r eference solar abundance was that of Grevesse et al 
(119961) . Once the best set of atmospheric parameters was cho 
sen for each star (see Section I3.2l i. we used the new MARC^I 
model atmospheres with standard compositiorOB- We chose the 
closest available global model metallicity (taking into account 
a-enhancement) of the ai Cen sub-populations, as reported in 
Tabled 

For all species we computed a 3 cr-clipped average of abun- 
dances resulting from each line. For elements that had both neu- 
tral and ionized lines, we computed the weighted (on the number 
of lines) average of the two ionization stages to obtain [El/Fe]. 
We typically rejected lines that had EW>100 mA, where the 
Gaussian approximation could fail, or EW<15 mA, since the 
relative error was too high. For elements with few lines, and 
where we had to rely mostly on strong lines, we either performed 
spectral synthesis, or checked that the DAOSPEC measurements 
were not too underestimated by visually inspecting the spectrum 
and overlaying the DAOSPEC Gaussian fit on each strong line. 

A few element abundances (Mn, Cu, Ba, and La) were de- 
rived with the help of spectral synthesis, taking into account HFS 
of single lines whe n need ed (see also Section ITTT i. We used the 
MOOCO ( Sneden,[l97l package in combination with the "best 
MARCS models" above to find the best fitting spectrurrQ All 



^ The broadenin g t heory adopted in this routine is that of 
lAli & Grieml ( Il965li and lVidal et atl il973h . 



We adopt the ob served EW rather than the th eoretical one (logg/ - 
d y., ■ lMagainLll98"4l) . following the discussion bv lMucciare 
' http://marcs.astro.uu.se/ 
This means [a/Fe]=+0.4 for metal-poor stars of [Fe/H]<-1.0 and 
reaching [Q'/Fe]=0 at [Fe/H]=0, following schematically the typical 
halo-disk behaviour of the Milky Way field population. 
' ' http://verdi.as.utexas.edu/moog.html 

To quantify the uncertainty related to the use of two different codes, 
one for the EW analysis and the other for spectral synthesis, we re- 
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Fig. 5. Examples of spectral synthesis applied to single lines 
for star WFI 512115. Big dots are the observed spectrum, solid 
lines represent the best fitting synthetic spectra, computed with 
all parameters fixed to the values of Tables [3] and |4] except for 
the abundance of the element under consideration. Spectra com- 
puted with the line abundances altered by ±0. 1 dex - with re- 
spect to the best fitting spectrum - are shown as dotted lines. 



3 0.5 




3BSE 
Wavalength (A) 




4305 



4310 

Wavelength (A) 



4315 



Fig. 6. Examples of spectral synthesis applied to the CN (top 
panel) and CH (bottom panel) molecular bands. Solid black 
curves represent the observed spectra for star WFI 512115, thick 
grey curves represent the best fitting synthetic spectra, while thin 
grey lines represent synthetic spectra differing from the best fit 
by ±0.2 dex in the N (top panel) and C (bottom panel) abun- 
dances. 



atmopheric parameters and the abundances of other blended el- 
ements (when present) were kept fixed and only the abundance 
of the element of interest was changed until the residuals of the 
fit were minimized. C and N abundances were derived by spec- 
tral synthesis of the molecular CN and CH bands at 3880 A and 
4300 A, respectively. Examples of line fits with spectral synthe- 
sis are shown in Figures |5] and |6l 

More in detail, carbon and nitrogen abundances were mea- 
sured by minimizing the residuals between grids of synthetic 
spectra and the observed ones around the CN and CH bands. 
In particular, C abundances were derived by fitting the CH G- 
band between 4300 and 4340 A, including the band heads of the 
(0-0), (1-1) and (2-2) bands of the A^A-X^H CH transitions; 
N abundances were derived with the CN B2E+-X22+ (0,0) UV 
band at 3870-3890 A. In the measurement process, we kept the 
atmospheric parameters and all the atomic abundances fixed to 
the values of Tables [3] and H) and we used the molecular data 
from the Kurucz database although, as discussed in Section lTTl 
a correction of -0.3 dex to the logg/ values of C was necessary. 
For each star, we computed synthetic spectra around the G-band 
by changing the C abundance only, and once the best fit was 
found, we employed the found C abundance in the fit of the CN 
band, where we changed only the N abundance. The procedure 
was then iterated a few times, repeating the fit for the CH band 
with the newly found N abundance, until the variations in both 
C and N were well below 0. 1 dex. 

The abundance calculation results are reported in Table |4] 
and are discussed in Section|4l 

3.4. Abundance uncertainties 

We estimated the internal, random uncertainty caused by imper- 
fections in the EW measurement (or spectral synthesis process) 
and in the line atomic data, as cr/ y/n, for those elements that had 
at least two surviving lines after a 3 cr-clipping pass. For ele- 
ments that relied on one line only, including those analysed with 
spectral synthe sis, we derived a typical uncertainty of 0.05 dex 
by means of the lCavrel(ll988h approximated formula, and we in- 
dicated this typical uncertainty between parenthesis in Table |4] 

The uncertainty owing to the continuum normalization pro- 
cedure can be estimated from the average spread of the resid- 
ual spectrum of each star after lines removal, which is automat- 
ically computed by DAO SPEC. For our spec tra, this was ~1%. 
According to Figure 2 by IStetson & Pancinol (|2008), this prop- 
agates to an approximate uncertainty in the EW estimates of 
+3 mA, and finally in an uncertainty of the order of 0.05 dex 
in the derived abundances. 

Another factor that has a big impact on the abundance ratios 
is the cho ice of atmospheric parameters (Section 13.2b . As dis- 
cussed bv lCavrel et al.l(l2004h . Teff, log g and v, are not strictly 
independent parameters when determined with the method of 
Section 13.31 Therefore, the best way to estimate the impact of 
the parameters' choice on abundance ratios is to change the most 
influential parameter, T^ff, and to re-optimize the other param- 
eters that naturally re-adjust to accommodate the temperature 
change. The difference between abundances calculated with the 
"best model" and with the altered one is a robust estimate of the 
systematic uncertainty owing to the choice of stellar parameters. 



analysed star WFI 507633 with MOOG, using the same linelist, EW, 
atomic data, atmospheric parameters and MARCS models. We found 
[Fe/H]=-2.02±0.02 which is well compatible with the abundance ob- 
tained with the Spite code (Table|4ll. 
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Table 4. Abundance ratios with random uncertainties (see text). 



Element 


WFI 503358 


WFI 503951 


WFI 507109 


WFI 507633 


WFI 512115 


WFI 512938 




(dex) 


(dex) 


(dex) 


(dex) 


(dex) 


(dex) 


fPel/Hl 


-1.70±0.01 


-1.58+0.01 


-1.99+0.02 


-2.05+0.02 


-1.18+0.01 


-2 03+0 02 


FFpTT/HI 


-1 70+0 17 


-1 56+0 OQ 


-1.95+0.07 


-1.94+0.05 


_i ^4+0 05 


_1 99+0.10 


TFe/Hl 


-1 70+0 01 


-1 58+0 01 


_1 QS+o 02 


-2 04+0 02 


-1.19+0.01 


_9 03+0 02 




+0 1 3+CO 05i 


j_n 40+0 7Q 


-0 1 7+l'0 051 


+0 14+0 50 


+0 80+0 05 


_0 19+('0 05i 

— yj' 1 y^yyj.yjj j 


TRa/Fel 




+ 1 1 0+0 24 


-0 06+0 24 


-0 20+0 1 6 


+0 64+0 1 2 


_0 S3+0 23 




-0 09+0 1 5 


-0 (S4+0 1 5 


+0 07+0 1 s 


-0.17+0.15 


-0.94+0.15 


_0 19+0.15 


rra/Fei 


+0 40+0 03 


+0 44+n 03 


+0 ^7+0 05 


+0 33+0 04 


+0 43+0 OS 


+0 34+0 03 


rrn/Fei 


-0 03+0 1 Q 


+0 71+0 1? 


-0 1 4+0 74 


+0 08+0 06 


+0 37+0 17 


-0 06+0 1 7 

— yj.yjyj^Kj. xz. 


rrr/Fei 


+0 OR+O 73 


-0 08+0 06 


-0 04+0 7 1 


_0 1 S+O 10 


+0 08+0 04 


+0 03+0 17 


rrii/Fpi 

vxj 1 CJ 




-0 60+fO 05i 


^_0 Q^+m 05i 


<-0 47+fO 05i 


__0 70+fO 051 


<_0 48+ rO 051 


n a/Fel 


+0 74+0 Ofi 




+0 17+0 08 


-0 1 3+0 1 


+0 1 8+0 09 


-0 06+0 08 


[Me/Fel 


+0.46+0.11 


+0.41+0.02 


+0.35+0.06 


+0.45+0.08 


+0.4I+(0.05) 


+0.55+0.10 


[Mn/Fe] 


-0.78+0.13 


-0.81+0.08 


-0.60+0.11 


-0.41+0.03 


-0.61+0.03 


-0.54+0.11 


[N/Fe] 


-0.12+0.30 


+ 1.36+0.15 


+0.57+0.20 


+ 1.23+0.15 


+ 1.26+0.25 


-0.29+0.30 


[Na/Fe] 


+0.07+0.11 


+0.60+0.12 


+0.30+0.03 


-0.15+0.33 


+0.67+0.19 


-0.42+0.37 


[Ni/Fe] 


-0.10+0.03 


-0.03+0.04 


+0.06+0.06 


+0.03+0.09 


-0.02+0.04 


+0.13+0.08 


[Sc/Fe] 


-0.03+0.07 


+0.20+0.05 


+0.10+0.12 


-0.07+0.09 


-0.04+0.12 


+0.18+0.14 


[Si/Fe] 


+0.38+0.14 


+0.46+0.30 


+0.27+(0.05) 


+0.19+(0.05) 


+0.35+0.39 


+0.61+(0.05) 


[Sr/Fe] 


+0.38+0.04 


+0.49+0.05 


-0.13+0.11 


-0.12+0.09 


+0.25+0.12 


-0.10+0.09 


[Ti/Fe] 


+0.25+0.04 


+0.38+0.03 


+0.30+0.04 


+0.33+0.06 


+0.38+0.04 


+0.24+0.03 


[Y/Fe] 


+0.58+0.08 


+0.71+0.06 


+0.31+0.05 




+0.69+0.13 


+0.28+0.72 



Table 5. Uncertanties owing to the choice of stellar parameters. Table 6. Comparative analysis of NGC 6397. 



Element 


WFI 512115 
(T=5350) 


WFI 507633 
(T=5650) 


UVES 
Average 


[Fel/H] 


+0.07 


+0.05 


+0.06 


[Fell/H] 


+0.05 


+0.08 


+0.07 


[Al/Fe] 


+0.02 


+0.02 


+0.02 


[Ba/Fe] 


+0.04 


+0.06 


+0.05 


[Ca/Fe] 


+0.03 


+0.04 


+0.04 


[C/Fe] 


+0.15 


+0.15 


+0.15 


[Co/Fe] 


+0.03 


+0.02 


+0.03 


[Cr/Fe] 


+0.17 


+0.04 


+0.11 


[La/Fe] 


+0.05 


+0.05 


+0.05 


[Mg/Fe] 


+0.05 


+0.02 


+0.04 


[Mn/Fe] 


+0.03 


+0.02 


+0.03 


[N/Fe] 


+0.20 


+0.20 


+0.20 


[Na/Fe] 


+0.01 


+0.02 


+0.02 


[Ni/Fe] 


+0.02 


+0.01 


+0.02 


[Sc/Fe] 


+0.06 


+0.06 


+0.06 


[Si/Fe] 


+0.02 


+0.02 


+0.02 


[Sr/Fe] 


+0.09 


+0.15 


+0.12 


[Ti/Fe] 


+0.02 


+0.05 


+0.04 


[Y/Fe] 


+0.05 


+0.06 


+0.06 



We therefore chose our warmest and coolest stars and re- 
computed their abundances with models having Teff altered by 
+ 100 K, re-optimizing the other parameters according to the 
method described in Section [331 The final uncertainties are ob- 
tained by averaging the absolute abundance differences of the 
-1-100 and -100 altered models, and are reported in Table |5] 

The global uncertainty (shown in all Figures fi'om[8]to[T3ll is 
computed as the sum in quadrature of the random eiTors, the un- 
certainty owing to the continuum placement and the one owing 
to the choice of atmospheric parameters. 



Star 






Vf 


A[Fe/H]'' 


Notes 




(K) 


(dex) 


(km/s) 


(dex) 




699 


5500 


3.4 


1.3 


+0.00 


Gratton et al. (20011 




5400 


3.3 


1.8 


-0.10 


Korn et al. (2007) 


793 


5500 


3.4 


1.3 


-0.05 


Gratton et al. (20011 




5500 


3.4 


1.6 


+0.06 


Korn et al. (2007) 


206810 


5500 


3.4 


1.3 


+0.09 


Gratton et al. (20011 




5500 


3.4 


1.8 


+0.07 


Korn et al. (2007) 


"Differences in the sense: ours minus literature, with same parameters. 



3.5. A comparison with NGC 6397 

Usually, to test the goodness of an abundance analysis, a compar- 
ison with the Sun or Arcturus is performed. In our case, because 
our linelist was optimized for metal-poor subgiants, we would 
be performing an abundance analysis of the Sun relying only on 
strong lines, with EW>100 mA. These strong lines are rejected 
when analysing our UVES metal-poor stars, because at this reso- 
lution they start to significantly deviate fro m the Gaussian shape 
(see Figure 7 bv lStetson & Pancinol l2008'). 

Therefore, while the general m ethod was already tested 
on the Sun by iPancino et al.l (l2010l) . we preferred to compare 
our present results to another well studied cluster, such as 
NGC 6397. We downloaded UVES archive spectra of three sub- 
giants (stars 206810 , 669, and 793) in NGC 6397, publi s hed by 
iGratton et aH (l200lh and later re-analysed bv lKorn et aP (l2007h . 
which have similar temperatures, gravities, and metallicities as 
the three most metal-poor SGB stars in our sa mple. We re- 
analys ed them adopting the sam e parameters as iGratton et al.l 
(I2001h and as iKorn et al.l ( l2007h . respectively, as detailed in 
Table |6l Obviously, the differences in the resulting [Fe/H] are 
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Rest Wavelength (A) 

Fig. 7. Comparison between the spectra of our UVES star 
WFI 507633 (topmost spectrum in the top panel) and star 
206810 in NGC 6397 from the UVES archive (bottom spectrum 
in the top panel). The bottom panel shows the difference between 
the two spectra. 



negligible within the quoted uncertai nties. We note in pa ssing 
that the major di fference between the iGratton et al.l (1200 ll) and 
iKornet al.l(r2007h abundance determinations for these three stars 
lies in the v, determination, wh ich is approximately 0.5 km s ' 
higher in the iKorn et al.l ( l2007h paper. This difference alone is 
probably enough to justify the different abundances found in the 
two studies. 

Figure Q compares the spectra of WFI 507633 and star 
206810 in NGC 6397. The parameters of these stars are very 
similar, and indeed the calcium and iron lines are virtually iden- 
tical, as confirmed by subtracting one spectrum from the other, 
where only noise remains. We are then co nfident that our analy - 
sis is correct and consistent with the one bv lGratton et al.l(l200lh . 
within the estimated uncertainties. 



4. Abundance results 

In the following sections we present and discuss abundance 
ratios of the measured elements. To compare with homoge- 
neous literature sets, we considered preferentially large (>10 
stars) samples of high-resolution measurements of stars be- 
longing to o) Cen. Unfortunately, there are no high-resolution 
studies of subgiantJp[ so we resorted to red giants surveys 
dNorris & Da Costal [l99l ISmith et al.L l2000t Ijohnson et all 
I2009h . This will help us to cross-identify sub-populations in 
the SGB region with those already identified on the RGB. 
For some species, there were no large literature surveys avail- 
able, so we included articles based on smaller samples (<10 
stars). The largest sample a vailable to date is the one by 
iJohnson & Pilachowskil (1201 Ol) . which contains more than 800 
stars, but at a moderate resolution (R^18 000). Their abundance 
ratio plots are more populous, but also obviously more scat- 
tered compared to higher-resolution studies, so we decided not 
to plot them in our Figures |8] to [13] except for their stars at 
[Fe/H]<-1.9 dex, which are not as well sampled in other less 
populous studies. We will also mention the abundace ratios by 



Although lVillanova et al.l ( 120071) study stars in the same evolution- 
ary phase as our targets, their resolution is lower (R^6000). 



IVillanova et al.l (l2010l) . a study of 38 stars observed in the exter- 
nal regions of cj Cen. Given the diversity of comparison studies, 
even if we reported all literature values to our solar abundance 
and \oggf system (using lines in common) when possible, resid- 
ual zeropoint differences among different studies will surely still 
be present. 

4.1. Iron 

Our abundance ratios for iron and iron-peak elements are shown 
in Figure [8] We note that three of our six stars lie around 
[Fe/H]^-2.0 dex, and that they do not correspond to any of 
the major sub-populations identified on the RGB, which are 
more metal-poor than any red giant studied with high-resolution 
spectroscopy before. The only exceptions are star ROA 213 by 
^mith et al. (2()( )0), w hich has [Fe/H]=-1.97 dex, star 85007 by 
fvillanova et all (l2010h at [Fe/H]=-1.98 dex, and a small group 
of 25 stars ( out of more than 800) below [Fe/H]^-2.0 dex in 
Figure 10 bv I Johnson & Pilachowskil (|2010|) . In their figure, the 
bulk populations have higher metallicities, with an abrupt drop 
in numbers below [Fe/H]^-1.9 dex, approximately. We tenta- 
tively classify our three stars as a separate subgroup - contain- 
ing few stars - that we name VMP (from very me tal-poor, see 
also S ection l5.1b . and we will consider ROA 2 1 3 b y Smith et all 
(l2000h and the 25 stars bv I Johnson & Pilachows ki (2010) as the 
VMP counterparts on the RGB in the following discussion. The 
other three targets apparently fall into the metal-poor and metal- 
intermedia te populations (MP and Mi nt, according to the classi- 
fication bv lPancinoetal .. 2000; Solli ma et all l2005al) . We will 
discuss population identifications in more detail in Section l572l 
We note that all comparisons in Figure [8] are between stars 
in different evolutionary phases (our targets are subgiants, while 
those in the literature are always giants). Therefore, the uncer- 
tainty on the zeropoint - especially as far as [Fe/H] is concerned 
- cannot be well quantified (see also Section l34l i. Possible fac- 
tors that cause systematically discrepant abundances between gi- 
ants and subgiants are 

- there is still no consensus for the NUTE corrections to 
Fel abundances, ranging from negligi ble ( Gratton et all 
119991) . to ap proximately 0.03-0 05 de x dKorn et all l2007l) 
to -1-0.3 dex dThevenin & IdiartL 1 1 9991) for stars of the type 
presented here; 

- the effect of diffusion could already be present for these sub- 
giants, lowering abundance by an amount that could be neg- 
ligible (IGratton et all 1200 Ih or of the order of 0.1-0.2 dex 
dKorn et al l 120071) ; however, diffusion is heavily influenced 
in all models by the amount of turbulent mixing, a phe- 
nomenon that is poorly constrained at present; 

- there are also possible granulation inhomogeneities, which 
apparently g o in the opposite direc tion compared to the 
above effects dShchukina et al.Ll2005h ; 

- finally, tw o studies that analysed these pheii omena in 
NGC 6397 dCiratton et al.ll2001l:lKorn etaill2007l) gave very 
different abundance results (see Section 13.51) . so even from 
the observational point of view it is quite difficult to provide 
reliable constraints to these phenomena. 

In spite of all these problems, we have the advantage that we 
do not use only [Fe/H], but we can also compare [El/Fe] abun- 
dance ratios - which should be more robust to measurement un- 
certainties and to variations caused by these physical processes - 
with a vast literature, to effectively cross-identify SGB and RGB 
populations. 
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Fig. 8. Abundance ratios of iron-peak elements. Our measure- 
ments are marked with filled black triangles: error bars in- 
clude random errors, continuum placement errors and errors 
owing to the choice of atmospheric parameters (see text). 
Grey symbols mark Uterature measureme nts for red giants: 
filled circles are from Johnson et ^ (120091) . empty pentagons 
from Cunha et al.l (1201 ol) , em pty circles fromlNorris & Da Costal 
( I995I). empty triangles fromlCohen' ("1981 ), filled squ ares from 



iGratton (Il982h. empty sta rs from iFrancois et al.l d 19881). asterisks 
from lCunha et al.l (l2002h . crosses from iJohnson & Pilachowskil 



(l20ldl) . four-pointed stars from lVillanova et alT i 20101) . 



4.2. Iron-peak elements 

While iNorris & Da Costal ( 1995h provided Cr a nd Ni measure- 
ments for 40 red giants, IJohnson et al] (l2009l) gave Ni mea- 
surements for 66 red giants and Cr for a smaller su bset. Cr 
and N i measurements are also pro vided bv |Villanova et al.l 
(I2OIOI) for 30-35 stars. We took the ICunha et al.l (l2002l) cop- 
per measurements for th e 40 red giants already analysed by 



iNorris & Da Costal (|1995|P^ For s t ars wi th [Fe/H]<-1.9 dex we 



used the Johnson & Pilachowskil (1201 Ol) Sc and Ni measure- 
ments. No literature measurements in large samples of red gi- 
an ts were found for C o and Mn, except for the recent Mn study 
bvlCunha et al.l(l2010l) . so we also compared with five giants by 
Cohenl(ll98ll). eight giants bv lGrattonI Jl982l) . and six giants by 



Francois et al1(ll988l) . 



There is general agreement for all iron-peak elements be- 
tween our measurements for SGB stars and the literature RGB 
ones. Two odd elements, Sc and Co, would give better results if 
hyperfine splitting (HFS) was taken into account, but since they 
do not give any additional information with respect to Fe, tra- 
ditionally no auhor has employed spectral synthesis to measure 
them. As a result, all ratios of Sc and Co in Figure [8]- including 
ours - are on average above solar Our Sc values fall exactly on 
top of literature measurements, and our Co values for the three 



'"^ T he Cu measurements of six stars by our g roup ([P ancino et"an 
are in substantial agreement with the ones bv lCunha et al.U2002l) . 
so we do not plot them in the various Figures. 




-2 -1.5 -1 
[Fe/H] 



Fig. 9. Abundance ratios of ff-elements. Symbols are the same 
as in Figure [8] 



most metal-rich SGB stars lie between the iGrattonI (IT982I) and 
the lFran9ois et alj (Il988l) measureme nts. It is known th at large 
NLTE corrections are required for Co dBergemanni |20 1 Ol) . which 
explains the rising trend observed in all cited studies. A slight 
underabundance can be noticed in the chromium ratios for the 
three most met al-poor stars, but N LTE effects should again be 
the cause for it dBergemanni |20 1 Ol) . It is interesting to note that 
[Ni/Fe] tends to be slightly higher for the VMP stars than for the 
remaining three UVES targets by 0.1- 0.2 dex, and that a simi- 
lar tre nd can be noticed in Figure 10 bv IJohnson & Pilachowskil 
d2010h . who do not comment on this effect. 

Mn and Cu are two very interesting elements with a dedi- 
cated vast literature because, although they belong to the iron- 
peak, they behave differently from other iron-peak elements (see 
iMcWilhaml [19971 for a classical review). In pa r ticular , Cu in 
w Cen has been first measured by iPancino et al 1 d2002h on SIX 
RGB stars, and later studied in detail bv ICunha et al.l d2002l) . 
who found it underabundant simila rly to field stars. Finally , 
it was theoretically modelled by Romano & Matteuccil d2007l) . 
who found massive stars as the m ost likely producers o f Cu, in 
agreement with previous studies dBisterzo et al.L l2004l) . While 
we find only upper limits for the three most met al-poor stars 

i Figur e |8]l, our measurements agree with those by ICunha et al.l 
2002h . 

Manganese instead was found to be extremely underabun- 
dant in our six stars. The determinations bv lCohenldl981h and 
Gratton ( 1982) appear only slightly subsolar, but none of these 
estimates takes into accou nt HFS. A bet t er agre ement was found 
with the measureme n ts by Cunha et al.l d2010l) . based on the 10 
stars by ISmith et akl d200C ), and derived with a complete HFS 
analysis, of which we plot the LTE resulting abundances in 
Figure |8] Our determinations, as those of Cunha etal. (20ldh. 
are lower than the values found for field stars dMcWilliaml 
1 19971) or metal-poor stars dCavrel et al.Ll2004 around [Fe/H]=- 
2.0 dex. In particular, three of the five analysed lines belong to 
the 40 00 A resonance triplet which, according to ICavrel et al.l 
d2004l) gives systematically lower abundances by approximately 
0.4 dex. Even rejecting these three lines, we found a [Mn/Fe] ^- 
0.7 dex for each of our six stars (Figure |8]l. The (marginal) dis- 
crepancy between our results and those by ICunha et alJ d2010h 
could arise because we used the Mn lines around 4000 A, while 
they use the lines around 6000 A. The general result - that 
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Fig. 10. Abundance ratios of heavy elements. Symbols are the 
same as in Figure [8] 



o) Cen has lower Mn abundance than field and GGC stars - 
suggests that the source of Mn production in aiC en should be 
low-m etallicity supernovae, either of type II or la (ICunha et all 

l2oToh . 



4.3. a-Elements 

It was not possible to measure the oxygen lines, or to put mean- 
ingful upper limits to the oxygen abundance in these subgiants. 
We were able to measure Mg, Si, Ca and Ti (Figure [9]l: all 
four elements were pre sent in the iNorris & Da Costal (fl995h 
and ISmith et alJ (l2000l) studies, but we plotted also t he Ca 
and Ti more recent measurements by iJohnson et al.l ( l2009l) 
and the Si, Ca, and Ti measureme n ts for stars with [Fe/H]<- 
1 . 9 dex by Johnson & Pilachowskil (|2010|) . The measurements 
bv lViUanovaet al.l (1201 Ol) also include all four a-elements stud- 
ies here. 



All measurements are slightly dispersed for Si because only 
a handful of lines is generally available, while for the other three 
elements they show a reasonably small spread. The three VMP 
stars do not seem to show any difference from the other three 
SGB stars in their ff-enhancement, with a weighted average of 
[a/Fe]=0.36±0.08 dex for all elements in all six stars. Literature 
measurements are quite scattered for Mg, but for RGB stars the 
Mg lines are very strong and difficult to measure, and Mg also 
anti-correrelates with Na, for example, so some spread should be 
expected. Apart from this exception, there is very good agree- 
ment among various literature determinations for giants, and 
with our measurements for subgiants. In particular, we note that 
RGB stars with [Fe/H]<-1.9 dex are always very close to our 
three VMP subgiants. The NLTE eff^ ects for all ou r SGB stars 
should be around ^0.1 dex in size dGratton et al.L Il999l) . for 
Mg, but not all the lines analysed here were also analysed by 
iGratton et alJd 19991) . so we prefer not to apply any correction. 

The (T-enhancement of VMP, metal-poor, and intermediate 
stars in a> Cen is typical of the field stars of the Galactic halo, 
which is indicative of a chemical enrichment dominated by 
SNe II. 



4.4. Heavy elements 

We were able to measure Ba, La, Sr, and Y. For Ba and 
La we used spectral synthesis (see Sections 13.11 and 13.31 1 
to take into account HFS. We plot the a bundance ratios in 
Figure [TOl along with t h e mea su rements by Norris & Da Costal 
(ll995h.nJ ohnson etalJ (l2009l) . IVillanova et alj ( |2010|) . and 
IJohnson & Pilach owskil (1201 (H) for stars with [Fe/H]<-1.9 dex. 



The log gf used by Norris & Da Costal d 19951) are now outdated, 
which is why, using the stars in common between literature stud- 
ies (i ncluding IVanture et all 11994 ISmith et all 120001: iPancinol 
l2003l) . we estimated that all their ratios needed to be raised by 
^0.4-0.5 dex - depend ing on the element - to c ompare with 
more recent studies. The lNorris & Da Costal ( Il995h data appear- 
ing in Figure [TOj are alr eady corrected f or this eff'ect. Also, the 
solar Ba abundance bv IVillanova et al.l (1201 Ol) was 2.31, while 
we use 2.13, and their data are corrected for this difference in 
Figure [To] as well. 

The heavy s-process elements Ba and La generally agree 
with literature estimates. In particular, the VMP stars ap- 
pear to have a lower average s-process enhancement, com- 
patible with zero, as confirmed by the La measurements by 
jjohnson & Pilachowski (2010). Thus, no s-process enrichment 
by AGB stars seems to have polluted the VMP star in the RGB, 
or in the SGB. This is supported by our Sr and Y measurements, 
which mainly agree with the literature measurements and which 
also have a low enhancement, compatible with zero for our VMP 
stars. 

It is also interesting that the most metal-rich star in our sam- 
ple, WFI 512115, appears to have a slightly lower s-process en- 
hancement than other stars measured in the literature at similar 
metallicity. If confirmed by larger samples of metal-rich stars, 
this could point out that the s-process enrichment by AGB stars 
in cl) Centauri was not completely h omogeneous, as supported 
also by the small group of stars in the lNorris &Da Costal d 19951) 
sample, around [Fe/H] ^-1.7, which appears to have [Ba/Fe] 0.2- 
0.3 dex higher than ot her stars at sim i lar me tallicity. A larger 
spread was found by IVillanova et al.l (1201 Ol) in s-process el- 
ements than in other studies presented here, which they ex- 
plain by assuming that a bimodality in these elements might be 
present at low metallicity, supporting the mentioned effect in the 
iNorris & Da Costal d 19951) data. Our three stars always share the 
same s-process enhancement, but larger samples are of course 
needed to see if the supposed bimodality, or higher spread, ex- 
tends to -2.0 dex stars as well. Finally, there is a large fraction 
of star s with solar Ba and Y enhancement in the Villanova et all 
d2010l) sample at all metallicities, which are not present in any 
of the other studies of red giants in o) Cen, and it is not clear yet 
if this is because of an intrinsic feature of stars in the outskirts 
of the stellar system, or a spurious measurement effect. 

4.5. Anti-correlations 

Figure [TT| shows the CH and CN band regions around 3880 and 
4300 A, respectively, for the three most metal-rich UVES stars. 
Clearly, two of them (namely WFI 503951 and 512115) show 
a deep CN band and almost no signal in the CH band. On the 
opposite, star WFI 503358 shows almost no CN band, and the 
CH band appears slightly deeper than in the other two stars. This 
evidence is not so clear in the three most metal-poor UVES stars, 
since the diatomic CN molecular band is much shallower 

We also have atomic lines of other elements that are found 
to (anti-)correlate within Galactic globular clusters, namely Mg, 
Na, and Al, besides C and N. Their abundance ratios are 
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Fig. 11. Smoothed spectra of the three most metal-rich UVES 
targets around the CN and the CH blue bands. 
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Fig. 12. Abundance ratios of Al, Na, C, and N for the six UVES 
stars. Symbols are the same as in Figure[8] 



plotted in Figures |9] and [12] An initial o bservation is, look 
ingat the literature a b undanc e ratio plots (Norris & Da Costa , 



1991 



Johnson et al 1 120091: iJohnson & PilachowskiL 120 lOt 

Marino et al.L l2010l) that [Al/Fe] and [Na/Fe] appear clearly bi- 
modal at all metallicities (with the possible exception of the 
most metal-rich stars around [Fe/H] ^-0.6 dex), with two dis- 
tinct sequences around [Al/Fe] ^0 and [Al/Fe] ^1 dex in the case 
of aluminium, and the same bimodality is visible, although less 
clearly, in the C and N ratios. In our Figure [T2l the Na data are 
confused because the lVillanova et al] (1201 Ol) data have an oppo- 
site trend to all the other literature data plotted, with Na decreas- 
ing with metallicity instead of increasing with it. We found no 
explanation for this trend. Apart from this, our data follow the 
literature trends reasonably well, with a similarly large scatter 
for [N/Fe]. We did not apply NLTE corrections following a rea- 
soning similar to that of Mg in Section 14.31 but thev should be 
again around 0. 1 dex for these subgiants dGratton et aUl200ll) . 

Figure [T3] reports on the usual (anti-)correlation plots. A 
clear and bimodal Na-Al correlation is seen in the literature data, 
as well as the C-N anti-correlation. Our data compare well with 
the literature in both cases. Even the most difficult Mg-Al anti- 
correlation is clearly visible in the literature data, although not 
as clearly in our own data for the six UVES stars. Finally, we 
also report the C-Al anti-correlation in Figure [T3] 

In summary, two of the three most metal-rich stars show 
clear signs of CNO pollution, as well as a clear Na-Al 
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Fig. 13. Anti-correlation plots for the six UVES stars. Symbols 
are the same as in Figure |8] 



anti-correlation, and this agrees with all past literature data 
on the subject. Our data for the three VMP stars are in- 
stead not conclusive about the presence of (anti-)correlations 
in the VMP population as a whole. We have two stars 
(WFI 507109 and 512938) with primordial composition and one 
star (WFI 507633) that seems to have high [N/Fe], but solar 
[C/Fe]. Literature measurements seem to point towards absent 
or reduced (anti-)correlations in VMP stars. The 25 stars by 
IJohnson & Pilachowskil(l2010l) with [Fe/H] <- 1.9 dex tend to all 
have high sodium, and a much lower dispersion than the other 
>800 stars in their sample. But there are no Al, Mg, C, or N mea- 
surements in their study. On the other hand. lMarino et aL 
show that for their red giants with [Fe/H] <- 1.8 dex, the Na-O 
anti-correlation appear less extended than that for other metal- 
licity groups in co Cen, but we do not know how the situation 
stands for [Fe/H]<-1.9 dex only. 

While a complete discussion on the anti-correlations among 
light elements in co Centauri is beyond the scope of the present 
paper - for the moment we just want to cross-identify our SGB 
stars with the RGB populations - we confirm that whatever the 
cause of anti-correlations in GC, it must have been active on 
each sub-population of o) Cen, possibly excluding only the most 
metal-poor one, at [Fe/H]^-2.0 dex, and the most metal-rich one 
(wher e all stars above -1 .0 dex appear polluted, iMarino et al.L 
I20T0I) . The presence or absence of (anti-)correlations in some 
sub-populations of u) Cen is a powerful tool to understan d its 
chemical evolution (see^Pa nciii oi 120031: ICarretta et al.LlIoTol and 
the discussion in Section lSTTT l. 



5. Discussion 

While comparisons between s ubgiants and giants are not entirely 
free from problems (see e.g., iBonifacio et al.L l2009l and refer- 
ences therein), we measured abundance ratios of several species, 
in good agreement with past literature determination for RGB 
stars. We found that three of our targets seem to belong to a sepa- 
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rate population with typical [Fe/H]^-2.0 dex, [a/Fe]-+0.35 dex, 
and an 5-process enhancement compatible with zero. Of the 
three remaining stars, two are consisten t with the abundanc e 
ratios of the RGB-MP (nomenclature by iPancino et all l2000h . 
and the last star a ppears to belong to th e RGB-Mint population 
(nomenclature bv lPancino et al.Ll2000t) or. more specifi cally, to 
the RGB-Mlnt2 (nomenclature bv lSollima et al.Ll2005ah . 

With this identification as our compass, we will try in the 
following sections to contribute to the understanding of the SGB 
of oj Cen on a few crucial topics. 

5.1. The first stellar generation in cj Cen? 



As discussed in Section 14.11 the three most metal-poor 
stars at [Fe/H]^-2.0 dex appear to belong to a separate 
sub-population, which we termed VMP (for very metal- 
poor). Stars as metal-poor as these were found in the 
past in small amounts, but were not considered as a sepa- 
rate sub -population per se. In particular, several (unbiassed 
with respect to metallicity) studies, aimed at deriving the 
metallicity distribution of oj Cen either through photometry 
jFrinchabov et al.L 120021: [SoUima et al.l l2005al: ICdamida et al 
20091) . low-resol ution sp e ctrosc opy of red giants ( Norris et al 



1996t ISuntzeff & Kraftt. [19961). or subgiants (iSolUma et al „ 
2005bl: [Stanford et al.L 120061: IVillanova et al.L 12007') found stars 



as poor as -2.0 dex or lower Some high-reso lution abun 
dance studies also fo u nd aga i n a few VMP stars (|Sm th et aT 



2000t iJohnson et al.L [20091: iJohnson & PilachowskiL I201C: 



Marino et al. , 2010t Villanova et al.L 20101) . and in general there 



appears to be an abrupt termination of the main MP population 
around [Fe/H]^-1.9 dex, with a sparse group of stars around 
[F e/H1^-2.0 dex. This effect is clearly visible in Figure 10 
M 



by [johnson & Pilachowskil (|2010() . This clean behaviour justi- 
fies the definition of VMP as a new, separate sub-population in 
(jj Cen. 

To estimate the fraction of such a minority population is not 
easy. From the cited studies we estimate that it should be at most 
5% of the entire stellar content of to Cen. Some more support 
to the existence of this small VMP component comes from the 
recent work by iBellini et alJ (1201 Ol) . who used exquisite ACS 
photometry to reveal additional sub-structure in the SGB re- 
gion of o) Cen. In part icular, the upper SGB branch (branch A 
in the nomenclature of 'Villanova et al. |j2007h appears split into 
two sub-branches in iBellini et al., (,20101) . The fact that three out 
of three of our upper SGB targets turned out to have [Fe/H]=- 
2.0 dex can thus be explained, because we chose them on the 
upper evelope of the upper branch, which turns out to be sepa- 
rated from the other branches, and which we can safely consider 
to be made of [Fe/H]=-2.0 dex stars, at least in the external re- 
gion of (o Cen that we are sampling here (see Figure|2]i. 

Our abundance ratios, based on high-resolution spec- 
troscopy, allow us to hypothesize that these stars must be the best 
candidate remnant of the primordial population in <jj Centauri, 
enriched primarily by type II SNe (given its [a/Fel^-nO.SS dex), 
and most probably free from severe pollution by AGB stars. 
This last statement is supported by the (almost) solar s- 
process ratios (Figure [TOli. Also, while three stars are too few 
to rule out the prese nce of anti-correlations in this popula- 
tion, literature studies (Smith et al. . 2000L Johnson et al l l2009t 



IJohnson & Pilachowskil2010t Marino et al. , 20101) suggest that 
(anti-)correlations should be absent or reduced in their exten- 
sion for VMP stars. It would be extremely interesting to study 
C, N, Mg, Na, and Al in larger samples of SGB stars in lo Cen. 
Indeed, although w Cen is commonly considered as the remnant 



of a dwarf galaxy accreted a long time ago by the Milky Way, all 
the sub-populations identified so far show clear anti-correlations 
(Figures[T2land[T3]) at all metallicities, except for - possibly - the 
VMP and RGB-a stars (iMarino et all 120101) . Field populations 
of dwarf galaxies do not normally show any anti-correlation , 
which is only found in globular clusters dGratton et al.L |2004|) . 
Therefore, the absence (or existence) of anti-correlations in the 
VMP component could give us evidence for (or against) the 
dwarf galaxy hypothes i s for the o rigin of u) Centauri, as dis- 
cussed also bv lPancind (l2003h and dCarretta et al.Ll2010l) . 

Summarizing, given the chemical properties of the VMP 
stars studies here and in the cited high-resolution studies, we 
conclude that they belong to a small and distinct sub-population, 
which appears to be the best-candidate (remnant) population of 
the first stellar generation in o) Cen, but could also - if the pres- 
ence of (anti)-correlations will be excluded - be the field popu- 
lation remnant of its hypothesized parent galaxy. 



5.2. SGB populations puzzle: who is who 

Apart from the notation used above, which separates the oj Cen 
sub-populations by metallicity in VMP, MP, Mint (in turn di- 
vided in Intl, I nt2, and Int3) , and MR (or RGB-a and S GB- 
a) proposed by iPancino et al.l (l2000h . iFerraro et al] (|2004 and 
SolUmaetal.l (20b5ai) . we will also use the nomenclature by 
Villanova et al. ( 20071) in this Section, who photometrically di- 
vide the SGB into four sub-branches (A, B, C and D) from the 
upper SGB envelope down to the SGB-a. This classification is 
also put into evidence in the two leftmost panels of Figure [T4| 

The whole difficulty in the study of the SGB of co Cen is to 
cross-correlate the five sub-populations defined photometrically 
and spectr oscopically a l ong th e RGB with the four or five (and 
more, see iBelUni et al.L 1201 Ot) sub-branches that are visible on 
the SGB. This is because, while on the RGB there are exquisite 
high-resolution studies that can link photometric features with 
chemical abundance patterns, this is not yet possible for the 
SGB. One important observation suggested by past attempts to 
link SGB and RGB populations is that the nicely combed struc- 
ture that appears in the RGB and the similarly nicely combed 
one that appears on the S GB could be somehow scrambled and 
mixed, as pointed out by IVillanova et al.l (|2007|) . While on the 
RGB the dominating factor that separates populations is metal- 
licity, on the SGB it must be age (together with C, N, O, and He), 
and this could complicate things significantly if there were no 
clear (monotonic increasing) age-metallicity relation in lj Cen. 
The uncertainties inv olved in the past low -res olution spectro- 
scopic studies such as lSoUima et al.l (l2005bl) and IVillanova et all 
(120071) were apparently not sufficient to give a final answer to 
the problem. The present study, on the other hand, while having 
higher precision in the abundance determination, can only rely 
on a limited number of stars. 



5.2.1. The VMP population 

Concerning our newly defined VMP population, we must as- 
sume that it should lie: ( i) on the bluest egde of the RGB colour 
distribution, because of its low metallicity, and (ii) on the up- 
per edge of the SGB, and i n parti cular on the upper edge of 
branch A bv IVillanova etal.1 (l2007h . where we find three stars 
out of three at [Fe/H]^-2 dex. This is supported by Figure 12 by 
Bellini et al. (2010), where for the first time the upper SGB ap- 
pears split into two separate branches, and where the lower base 
of the RGB shows a clearly bluer sub-sequence made of a small 
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Fig. 14. An example of isochrones fit (see text) to determine age differences among the three sub-populations studied in this paper. 
Left panel: the literature values of E(B-V) and (m-M)v for u Cen are used to derive age differences; the oldest populations turns 
out to be 16 Gr old. Centre panel: a cosmological age of 13 Gyr is forced for the oldest populations, altering E(B-V) and (m-M)^ 
but with no significant effect on the resulting age differences; the shape of the isochrones is now in worse agreement with the data. 
Right panel: NGC 6397 is shown for comparison. The letters A, B, and C, wh ich appear in the two leftmost panels, indicate the 
approximate location of the SGB branches as defined bv lVillanova et al.l (l2007h . 



number of stars. These features are the most likely photometric 
counterparts of the VMP population found here. 



5.2.2. The MP population 

The connection between the SGB-MP and the RGB-MP is not 
clear yet. While some other population should be present to- 
gether with the VMP to account for the number of stars in the 
upper branch, or branch A, it is not clear if this is really the 
counterpart of the RGB-MP, because the RGB-MP contains at 
least 30% (and likely more) of the stars in oj Cen. We cannot 
speculate too much with the data in hand - the only certain fact 
is that the two MP stars in this paper (WFI 503358 and 503951) 
both lie lower than branch A, possibly on branch B or even C - 
but we see two possible solutions: 



5.2.3. The Mint populations 

Our star WFI 512115, which chemically belongs to the Mint 
population (and, more specifically, to the MInt2) also appears 
to lie either on branch B or - most probably - on branch 

C. Combining this e vide nce with the one pre sented by both 
ISolhma et all i2005bl) and lVillanova et al.l (l2007l) , we tentatively 
conclude that the most likely photometric counterpart of the 
MInt2 population should be branch C. The presence of a few 
RGB-MP stars on branch C stiU cannot be ruled out, and it 
could be either caused by some scatt er in the low-resolution 
abundances by iVillanova et al.l (l2007l) . o r by photome t ric er- 
rors in the W FI reference photometry jPancino et al.L 120001; 
iPancinoi l2003h, or even by a different SGB shape for the SGB- 
MP (iBellini et all 1201 Oh . We also speculate that the RGB-MInt3 
should correspond to a sequence lying between branches C and 

D, as clearly indicated by .Soll im a et al . (2005b) in their Figure 4, 
and mentioned in passing bv lvillanova etal.l(l2007il . 



One possibility - supported by our analysis of stars 
WFI 503358 and 503351 - is that the RGB-MIntl stays 
on branch A (lower half) together with the VMP, and the 
RGB-MP stays instead on branch B. This would allow for a 
be tter agreement between the population fract i ons es timated 
bv ISoUima elaP (l2005ah and I Villanova et al.l (l2007h . while 
the uncertainties of both low resolution studies are probably 
large enough to accommodate such a switch. 
Th e other obvious interpre tation has been already mentioned 
by IVillanova et al.l (l2007l) . i.e., that MP stars He on both 
branches A (lower) and B, mixed with part of the metal- 
intermediate populations. 



5.2.4. The MR or SGB-a population 



While the 
(iPancino et al.l 



c ross-i dentification 
2000t 



between 
2005ah. 



the 
the 



RGB-a 
SGB-a 



iSoUima et al.L . 

'Ferraro e t al!r(l2004;^[SoUima et al.l (l2005bl) and branch D 
(Villano va et al. n2007h appears unquestionable, there is still 
some debate about its metallicity, ranging from [Fe/H]^-l.l, 
derived from low-res olution spectroscopy of SGB stars by 
IVillanova et"al] (|2007|) to [Fe/H]^-0.6 de x, derived from both 
high-resolution spectroscopy of RGB stars (Pancino et al.','200] 
and lo w-resolution spectroscopy of SGB stars (Sollima et al 
l2005bh . The metahicity of the SGB-a wih be the subject of a 
following paper based on GIRAFFE data, but it is interesting to 
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note here that lBellini et al.l (1201 Oh found the SGB-a (or branch 
D) split into two sequences, suggesting that two populations 
with slightly different properties could occupy these two very 
close loci. 

5.3. Towards a solution of the age spread problem 

The history of age spread determinations in w Cen comprises a 
variety of studies, all focused on the SGB, the most age sensitive 
region in the CMD. An age spread of 2- 5 Gyr, among the var- 
ious s ub-po pulations was first found by iHu^hes & WallersteinI 
(120001) and Hilker & Richtlej (l2000l) . based on the TO region 
morphology of their high-quality Stromgren photometries and 
colour-metallicty calibrations. Later several other papers came 
out based on combinations of high quality photometry and low- 
resolution spectroscopy (Hughes e t al., 2004; H ilker et al ., 2004 : 



RevetaLL l2004 LSo lima et al.L l2005bt LStanfordet al.[ l200 



Villanova et al .L l2007h . finding again various age dispersions all 



around 2-5 Gyr The only studies repor t ing on age difference 
below 2 Gyrs are those by iFerraro et al.l (l2004l) . who find that 
the SGB-a cannot be fitte d with any isochrone y ounger than the 
SGB-MP population, and lSolhma et al.l (l2005bl) . who similarly 
found that the overall age spread of the SGB sub-populations 
cannot amount to mo re th an 2 Gyr It is in t erestin g to note that 
ISoUima et al] (l2005bh and IVillanova etal] (|2007|) use photom- 
etry from the same ACS data set and s pectra of similar qual- 
ity, and reach opposite conclusions, with IVillanova et all ( |2007|) 
finding an age spread of at least 2 Gyr within the MP population 
only. 

On the one hand, spectroscopic abundances of both RGB and 
SGB stars find an a-enhancement of all populations (except the 
RGB-a) consistent with pure type II SNe enrichment which, in 
the self-enrichment scenario for co Cen, would imply very fast 
enrichmemt, within 1 Gyr On the other hand, the high .^-process 
enhancement of all populations (except maybe the VMP found 
here) would imply enrichment by intermediate mass (1-3 M©) 
AGB stars on longer ti mescales (1-2 Gyr and possibly even 
more lBusso et al.Lll999l) . Therefore, if oj Cen is a self-enriched 
system, we should find (and indeed many authors do find) some 
significant age spread. While we cannot give the ultimate so- 
lution to the relative ages (or age spread) puzzles listed above 
with the present data, we can nevertheless use our newly defined 
VMP population to shed some light on the problem. 

We use isochrones from the Padova database of stellar evo- 
lutionary tracks a nd isochronern, and in pa rticular we chose the 
ones based on the lMarigo et aljT2008h andlGirard i et al.l (I200C ) 



tracks . We preferred this set over BaSTO (Pietr inferni et al 
|2006|) simply because they provide tranformations to the WFI 
filters. We note that the results in ages do not change signif- 
icantly when using BaSTI (less than l^Gyr overall shift), but 
the actual shape of the SGB is better reproduced when the WFI 
filters are used instead. While the absolute ages are somewhat 
uncertain, we show in Figure [14] that the age differences are rel- 
atively robust regardless of the absolute age scale calibration of 
the chosen isochones set. 

Firstly, we computed the Z values for the three populations 
using the formula 

logZ = [Fe/H] - 1.7 -H log(0.638 10 ^"'^""^ + 0.362) 

according to ISalaris et al] (Il993h. We adopted the iLubl (l2002h 
reddening and 



Bellazzini et al. (|2004|) distance modulus (see 



also Section 13.21 ). We used the cosmological helium abun- 
dance models (for a discussion on higher He abundances see 
Section |54] |. The uncertainties involved in the isochrone fitting 
procedures are large: we assume aproximately ±2 Gyr. We find 
that the VMP and MP populations must be relatively close in 
age, with a difference of 0+2 Gyr, while the Mint component 
should be about 2+2 Gyr younger. If we force the age of the 
oldest population(s) to a cosmological value of 13 Gyr, we have 
to increase E(B-V) and {m-M)v by 0.03 and 0.16 mag respec- 
tively, and the fit becomes considerably worse when consider- 
ing the isochrones shape. The age differences, however, do not 
change. As a comparison, we show in Figure [14] the case of 
NGC 639 7, using the WFI B, V and I photometric catalogue de- 
scribed bv lCarretta et alj 
modulus from the revised 



20091) and the reddening and distance 
Harris] (1 1996h catalogue. 



http://pleiacli.pd.astro.it/ 
www.oa-teramo.inaf.it/BASTI/ 



The existence of the VMP at the SGB level could help in 
solving two of the paradoxes previously found in the literature. 
The first is the zero age differe nce - within the uncertainties - 
found by ISoUima et all (l2005bl) between the MP and the Mint 
populations, which is difficult to understand given the difference 
in s-process enhancement between the two populations. If one 
assigns [Fe/H] =-2 dex to the upper SGB envelope, and moves 
the MP population to branch B, this shift in metallicity wou ld be 
allowed by the uncertainties in the Sollima et al.' ('2005b) low- 
resolution spectroscopy (see also Starkenbur g et al., 2010, on the 
uncertainties of the calcium triplet calibration at low metallicity). 
As a result, the MP would become older, and the age difference 
between the MP and Mint populations would not be so close to 
zero anymore. 

The second puzzle that could be alleviated by the existence 
of the VMP is the coexistence within the MP group of two sep- 
arate populations with ages differing by ^2-3 Gyr, found by 
]Villanova et al.] (]2007]) . which could be explained if one admits 
that the MP population is in reality a mix of MP and VMP 
stars. The metallicity difference between the VMP and MP pop- 
ulations would be difficult to resolve with low-resolution spec- 
troscopy, si milarly to the ca s e abov e. The younger MP group 
identified by ]Villanova et al.] (]2007h would become older if its 
metallicity was [Fe/H]^-2.0 dex instead of -1.7 dex, thus recre- 
ating an almost monotonic age-metallicity relation, although 
the age difference between the VMP and the MP populations 
would be smal l . It co uld be argued that the younger MP in 
]Villanova et al.] (]2007l) is much more numerous than the VMP 
estimated fraction, i.e., 5% of the total stellar content in co Cen. 
However, we must note that a large fraction of the targets in 
that study are selected a priori on the upper SGB branch (or 
branch A), which we suspect is dominated by VMP stars. In this 
case, the target selection would be biassed preferentially towards 
the VMP popu l ation and the relative numbers in Figure 19 by 
]Villanova et al.] (]2007]) would not be representative of the respec- 
tive population fractions anymore. 

In a self-enrichment scenario, the small age difference be- 
tween the VMP and the MP population (0+2 Gyr) poses no prob- 
lem as far as the type II SNe are concerned, but could perhaps 
be too short to accommodate the ^0.5 dex overabundance in s- 
process elements for the MP. On the other hand, if 2±2 Gyr oc- 
curred between the MP and the MInt2 population, there could be 
enough time to enrich the Mint population in s-process elements 
up to the observed level of [s/Fe]^l dex. Detailed chemical evo- 
lution calculations would be extremely useful in understanding 
these details. 
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5.4. The helium abundance 



Since the di s cover y of a double MS in to Cen (lAndersonll2002l: 
iBedin et all |2004|) . and the evidence that the bluer sequenc e 
was more metal-rich than the redder one dPiotto et all l2005h . 
it was suggested that some of the metal-ric her population s in 



a> Cen could have abnormal He abundan ce dNorris 



2001 . of 

about Y^O.35-0.40. A recent review by iRenzinil (l2008l) dis- 
cussed the possible scenarios and compared them with obser- 
vations of u) Cen, b ut also of NGC 28 08, which was found to 
possess a triple MS (iPiotto et all l2007h . along with other mas- 
sive GGC. 

Such an overabu ndant helium would have some impact in the 
model atmospheres dBohm-Vitensei 119791: iGirardi et all l2007h 
that are at the basis of any abundance analysis such as the one 
presented here. The first thing to note is that the most metal-poor 
populations of to Cen do not require any He enhancement, and 
we assume this to be the case not only for the two MP stars 
(WFl 503358 and 503951), but also for the three VMP ones 
(WFl 507109, 507633, and 512938). The only star that could 
suffer from an helium enriched atmoshpere is the one belonging 
to the Mlnt2 population, WF l 5121 15; 

A simplified treatment by |Gray| ( l2008h assumes that an over- 
abundance in helium has a similar effect as an increase in gravity, 
as a first approximation: 

^g 4 AA(He) 



g 1 -H 4 A(He) 

Just to derive an order-of-magnitude effect, we translate a 
mass fraction increase from Y^O.25 to Y^O.35 into a number 
increase from A(He)^0. 10 to A(He)^0. 15. This would be mim- 
icked, in our WFl 512115 star, by an increase in surface grav- 
ity from roughly 3.5 to 4.0 dex, corresponding to an increase 
in abundance of 0.08 dex in [Fe/H]. Thus, WFl 512115 would 
change from [Fe/H]=-1.19 to -1.11 dex, with basically no sig- 
nificant effect in the age difference determination, within the 
quoted uncertainties. We finally note here that, as discussed by 
ISoUima et al.l (l2005bl) in their Figure 6, this increase in the he- 
lium abundance would only change the shape of the isochrone's 
SGB, making it steeper, with a negligible (less than 1 Gyr) im- 
pact on the relative age determination. 



6. Summary and Conclusions 

We analysed UVES high-resolution spectra of six stars on 
the SGB of u) Centauri. We compared our results with RGB 
high-resolution spectroscopy to identify the sub-populations to 
which our targets belong, and we found remarkable agree- 
ment with past abundance determinations. Three of our tar- 
gets (WFl 507109, 507633, and 512939) have [Fe/H]^-2.0 dex, 
are (^-enhanced and show no significant s-process enhancement. 
Two of the remaining targets (WFl 503358 and 503951) belong 
to the MP population, with [Fe/H] ^-1.65 dex, c-enhanced and 
with [s/Fe]^-i-0.5 dex. The last target, star WFl 5121 15, belongs 
to the Mlnt2 population, with [Fe/H] =-1.19 dex, a-ehanced 
and with s-process enhancement similar to the MP targets, i.e., 
slightly lower than what is expected for Mint stars. 

Our main result (see Section 15.11 ) is that there exists an 
additional, metal-poor population (that we name VMP) at 
[Fe/H] ^-2.0 dex, which has chemical properties that make it 
the ideal candidate (remnant of) th e primordia l popu lation 
of u Cen The RGB sta r ROA 213 ISmith et al.l dioOOl) . star 
85007 by IVillanova et"an d2010l) . and 25 red giants studied by 



iJohnson & Pilachowskil d2010h . have similar chemical composi- 
tion and could represent the prototype VMP members along the 
RGB. In particular, the s-process enhancement of the SGB- VMP 
population is not compatible with the RGB-MP stars, while it 
is compatible with the quoted RGB-VMP stars. Our conclu- 
sion is also supported by previous work on metallicity distri- 
butions or R GB and SGB stars and by the exquisite photom- 
etry by Belhni et all d201(]h . We estimate that this VMP popu- 
lation should comprise at most 5% of the entire stellar content 
of o) Centauri, at present. The presence or absence of light ele- 
ment anti-correlations in this population would be a fundamen- 
tal constraint to the nature of to Cen, because anti-correlations 
are generally exclusively found in globular clusters and never 
in the field populations of galaxies. From the available litera- 
ture (mainly | Johnson et alll2009l : |johnson & Pilachowskill20ldl : 



iMarino et all 1201 Ol) it appears that (anti-)correlations could be 
reduced in extent, in VMP stars. Until the presence of (anti)- 
correlations in VMP stars is excluded by larger data samples, 
it looks more promising to interpret this as the primordial pop- 
ulation of (jj Cen instead of the remnant field population of its 
putative parent galaxy. 

The high-precision abundance determinations obtained al- 
lowed us to try to shed some light on the SGB morphology rela- 
tion with the RGB spectroscopically identified sub-populations. 
We conclude that there appears to be no one-to-one correspon- 
dence between the nicely combed substructures of the RGB and 
SGB. In particular, the MP and Mint populations could either 
be mixed along the (lower) A, B, and C branches, or be posi- 
tioned in a not strictly monotonic order, with Mlntl occupying 
branch A (lower) and MP branch B, just as an example. As al- 
ready said, VMP stars should occupy (and pos sibly dominate) 
the up permost SGB branch, or branch A in the IVillanova et al.l 
d2007h terminology. 

We also found that (see Section l53] ) the existence of the VMP 
population could alleviate some of the problems found in pre- 
vious determinations of relative ages. In particular, the small 
metallicity difference between the VMP and MP populations 
could have escaped previous abundance analyses based on low- 
resolution spectra. The puzzling result byfVillanova et al. (2003) 
that the MP population should contain two groups with different 
ages could indeed be explained by the metallicity difference be- 
tw een VMP and MP. Al so, the (too) small age difference found 
bv lSoUima et alj d2005bl) between the MP and Mint populations 
woul become slightly larger when taking into account the exis- 
tence of the VMP, which should dominate the uppermost SGB 
envelope. 

Finally, there should be a small age difference between the 
VMP and MP populations ( 0+2 Gyr), while a slightly larger age 
difference (2+2 Gyr) should occur between the VMP and the 
MInt2 populations. Althought this latter result is less secure be- 
cause it relies on one star only , it agrees very weU with the ma- 
jority of past studies (see, e.g.. [Stanford et aLLl2006h . The use of 
different sets of isochrones (Padova, BaSTI) does not change the 
result significantly, and the helium abundance problem should 
have a negligible impact on the Mlnt2 star WFl 512115 (the 
only one which should have higher helium) because, to first ap- 
proximation, it should change its [Fe/H] by ^0.08 dex. The age 
distribution suggested by the present data would accommodate 
a fast enrichment between the VMP and MP populations, domi- 
nated by SNe type 11, while the s-process enrichment of the MP 
([s/Fe]^+0.5 dex) could still pose a problem. The age difference 
between the MP and Mint populations could instead be suffi- 
cient to allow for so me intermediate-mass AGB star enrichment 
dBusso et all 1 19991) . bringing [s/Fe] to +1.0 dex. 
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We conclude by noting that this is the only high-resolution- 
based abundance analysis published on SGB stars in co Cen so 
far Even if the precision of the abundances is higher than in past 
low-resolution studies, of course the number of stars examined 
is only six. To give the final answer to the relative ages problem 
in CO Cen, and to identify who is who in the CMD at the SGB 
level, a much larger sample (a few hundreds) of relatively high- 
resolution spectra in the SGB region is absolutely necessary. 
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